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Abstract. In a statistical analysis of Debrecen Photoheliographic Results sunspot area 

data we find that the logarithmic deviation (logD)' of the area decay rate D from the 

parabolic mean decay law (derived in the first paper in this series) follows a Gaussian 

probability distribution. As a consequence, the actual decay rate D and the time-averaged 

decay rate D are also characterized by approximately lognormal distributions, as found in 

qq ■ an earlier work. The correlation time of (log D)' is about 3 days. We find a significant phys- 

iy\ | ical anticorrelation between (logD)' and the amount of plage magnetic flux of the same 

£SJ ■ polarity in an annulus around the spot on Kitt Peak magnetograms. The anticorrelation 

\Q ' is interpreted in terms of a generalization of the turbulent erosion model of sunspot decay 

to the case when the flux tube is embedded in a preexisting homogeneous "plage" field. 

The decay rate is found to depend inversely on the value of this plage field, the relation 

being very close to logarithmic, i.e. the plage field acts as multiplicative noise in the decay 

process. A Gaussian probability distribution of the field strength in the surrounding plage 

will then naturally lead to a lognormal distribution of the decay rates, as observed. It is 

thus suggested that, beside other multiplicative noise sources, the environmental effect of 

surrounding plage fields is a major factor in the origin of lognormally distributed large 

random deviations from the mean law in the sunspot decay rates. 
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1. Introduction 

The very large and apparently random individual differences in the area 
decay rates of sunspots have constituted a serious obstacle in the way of 
an empirical determination of a mean decay law. With time it has been 
recognized that the study of the distribution of decay rates itself may be 
of key importance for the proper statistical analysis of the decay, as most 
of the widely applied tools of mathematical statistics (e.g. least-square fits) 
are based on the tacit assumption of one particular (usually the Gaussian) 
probability distribution. The realization that the mean area decay rates D of 
sunspots are distributed lognormally (Martinez Pillet, Moreno-Insertis, and 
Vazquez, 1993, hereafter MMV93) was therefore an important step towards 
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the proper understanding of the statistical regularities underlying the ap- 
parently haphazard sunspot decay process. A lognormal distribution implies 
that it is the logarithm of the decay rate D that is normally distributed, and 
thus it is log D that should be used as the dependent variable in least-square 
fits and other standard statistical procedures. This realization has helped 
us in the first paper of this series (Petrovay and van Driel-Gesztelyi, 1997, 
hereafter Paper I; see also Petrovay, 1998) to finally determine, by rigorous 
statistical methods and at a quite convincing confidence level, the mean law 
governing sunspot decay. It was found that an "idealized" sunspot following 
this mean law exactly, with no random deviations, would decay according 
to the law 

Adcal = C D r/r C D = 32.0 ± 0.26. (1) 

where -Didcal is the area decay rate D for an idealized spot in units of 
MSH/day (MSH: one millionth solar hemisphere), r is the equivalent radius 
of the spot (i.e. the radius of the circle with the same area as that of the 
spot), while ro is its maximal equivalent radius. (The corresponding areas 
will be denoted by A and Aq, respectively.) It is clear that as D = —A and 
r ~ A 1 ' 2 , Equation (1) implies a parabolic decrease of the spot area with 
time, in contrast to the previously widely held view that the decay is linear. 
For real spots the actual daily decay rate D shows random deviations 
from the law (1). It follows from the considerations of the first paragraph 
above that in the derivation of the law (1) it is the logarithm of D that had 
to be averaged, i.e. the mean decay law for real sunspots is 

(logD>=log(C D r/r ) (2) 

(Here we are using the notation (a) for the ensemble average over many 
spots at the same phase of decay r/ro, as opposed to the time average over 
the decay phase of one spot, a. Deviations from the ensemble average are 
denoted as a' = a — (a).) 

Having determined the mean decay law, in the present paper we turn our 
attention to the deviations from this mean law. There are several important 
issues to clarify in this respect. 

- The analysis of Paper I assumed that (log D)' follows a Gaussian dis- 
tribution. In the case of a parabolic decay law like (1) however it is 
not possible to have a strictly normal distribution simultaneously for 
(log-D)', and also for logD and logD, as suggested in MMV93. In 
Section 2 we therefore consider the question, which one of the pro- 
posed lognormal distributions is the fundamental and exact one, and 
which are those that are only approximate and appear as secondary 
consequences of the fundamental distribution. 

- The decay law (1) implies that the decay process has a "memory": 
in order to "know" the "right" decay rate, the spot must remember 



decay2pp.tex; 1/02/2008; 10:30; p. 2 



Making Sense of Sunspot Decay, II o 

its maximal radius r^ at all later phases of its decay. The question 
arises, how can this long-term memory be reconciled with the presence 
of random fluctuations in the decay rate? In Section 3 we first estimate 
the correlation time of the fluctuations, and then, considering the pos- 
sible explanations, we conclude that the "memory" responsible for both 
the long-term systematics and for the short-term fluctuations can be 
identified with the amount of flux stored in the plage around the spot. 
This implies that decay rate fluctuations should be correlated with the 
amount of this flux. 

- In order to check the above prediction, in Section 4 we compare the 
amount of magnetic flux in an annulus surrounding the spot, as de- 
termined from Kitt Peak magnetograms, with the decay rate. We 
indeed find a significant anticorrelation between the two quantities. 
Interpreting this empirical anticorrelation, in Section 5 we generalize 
the turbulent erosion model to account for a plage field. 
Finally, in the Conclusion we discuss and summarize the main factors 
leading to, and the physical interpretation of the fluctuation in the decay 
rates. 

Throughout the paper, we use the same data set as in Paper I, taken 
from the Debrecen Photoheliographic Results (DPR) 1977-78 (Dezso, Ger- 
lei, and Kovacs, 1987, 1996). This consisted of 3990 area measurements 
of 476 different spots. For a detailed discussion of the data selection and 
reduction process we refer the reader to Section 3 of Paper I. In order to 
extend the data set, we made an attempt to include data from the more 
extensive Mt Wilson catalogue. However, owing to the lower precision of 
those data and the problems with the day-to-day identification of sunspots 
we decided to use the DPR data only for the present analysis. Nevertheless, 
the comparison of the two catalogues offered some interesting conclusions, 
summarized in Appendix A. 



2. The distribution of decay rate fluctuations 

Figure la shows the histogram of (log-D)', as defined after Equation (2), 
for our sample (dash-dotted line). This original histogram is distorted by 
the selection effect described in Section 3.3.2 of Paper I: this effect essen- 
tially consists of a lower expected number of observations for spots with 
higher decay rates, and thus shorter lifetimes. (The histogram contains 
data points with r/ro > 0.7 only, to avoid the regime most affected by 
this bias.) Correcting for the bias using the method described there, we 
obtain the histogram shown by the solid line. A Gaussian fit to this his- 
togram (dashed) is found to be a fairly good representation of the data. 
(Q = T[(N - l)/2, x 2 / 2 ] = 0-02; tne skewness and kurtosis are -0.20 ± 0.21 
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Figure 1. Histograms of (a) (logD)' and (b) logD for all spots that died on the visible 
hemisphere. (Fig. la contains data for r/ro > 0.7 only.) Dash-dotted: raw histograms; 
solid: histograms corrected for selection effect; dashed: Gaussian fits. The parameters of 
the fits are ((log-D)') = -0.09, a = 0.19 for (a) and (logD) = 1.25, a = 0.23 for (b). 




Figure 2. Sketch for the explanation of the relation between the lognormal distributions 
of different decay rates. S is a generic data point in the plane r/ro-log D where the thick 
curve represents the mean decay law (1). The thin curve is a decay curve for a generic 
sunspot with a random logarithmic deviation (logD)' varying on a timescale r, while 
the dotted curve is the decay curve of a hypothetical suspot with a constant logarithmic 
deviation from the mean law, (logD)' =const. 



and 0.67 ± 0.52, respectively, i.e. not significantly different from 0.) This 
normal distribution of (logD)' confirms the correctness of the analysis of 
Paper I, based on least-square fits with log D as the dependent variable. 

In MMV93 it was found that, within the errors, \ogD is normally dis- 
tributed. Furthermore, based on a statistics of Mt Wilson data (Howard, 
1992) it was suggested that log D should also show a Gaussian distribution. 
However, our above finding on the normal distribution of (log-D)', together 
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with the parabolic mean decay law (2), implies that the distributions studied 
in MMV93 cannot be exactly Gaussian. 

To see this first for log.D = (log-D) + (log.D)', let us consider the plane 
r/rQ-logD (Figure 2). The normal distribution of (log-D)' with zero mean 
implies that, for a narrow range of r/ro (shaded area), log-D is also normally 
distributed around a mean value log(Cx>r/ro). This mean depends on r/ro; 
hence, the distribution of log.D in the complete sample (essentially, the 
distribution of all data points in the plane projected on the ordinate) will 
deviate from a Gaussian. (This is in contrast to the case of linear decay 
where the mean decay law, indicated by the thick curve on Figure 2, would 
be a horizontal line, independent of r/ro.) Indeed, in Figure lb we present 
the histogram of log-D: while at first glance the quality of the fit does not 
seem to be very different from the previous case, a quantitative analysis 
shows that the Gaussian fit is insufficient to explain the data. (Q < 10~ 5 ; 
the skewness and kurtosis are — 0.77±0.13 and — 0.13±0.33, respectively, i.e. 
the skewness significantly differs from 0.) The deviation from gaussianity is 
relatively mild owing to the fact that the number of data points is higher for 
higher values of r/ro (cf. Fig. 1 of Paper I) where the thick curve representing 
the mean law is not far from horizontal. The distribution of log D is therefore 
indeed approximately lognormal, as proposed in MMV93; but, as we pointed 
out, the deviation is statistically significant. 

Let us now turn our attention to the time-averaged decay rate, defined 

as 

— 1 f t2 

D = - -/ Ddt. (3) 

*2 - *1 Jt! 

Defining logC^, = logCo + (l°g-D) / ) according to Equation (2) we may 
formally write 

D = C' D r/r . (4) 

For a hypotetical sunspot with (log-D)' =const. (dotted curve in Figure 2) 
C' D =const., so from Equation (4) 



log D = log C' D + log r/r (5) 

where the last term is a constant for all spots, as for a parabolic decay law, 
r decreases linearly with time, (r/ro = 0.5 if t\ and £2 correspond to the 
area maximum and the total disappearance of the spot, respectively.) So for 
such spots if (log-D)' is normally distributed, so is logC^ and log-D. For 
real sunspots, however, (log-D)'=const. does not apply; instead, their actual 
decay rate will fluctuate around the mean law with some timescale r (thin 
curve in Figure 2), and the above reasoning does not hold. Note that if the 
mean decay law were linear ((log-D) =const.), the resulting distribution of 
log -D would still be normal, albeit with a lower dispersion a, as in that case 
the last term in Equation (5) would not be present. 
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Thus, D in general cannot be exactly lognormally distributed; but its 
actual distribution may still appear to be close to lognormal provided either 
that (log D)' does not change very much in the time interval [£1, £2] (i- e - that 
£2 — £1 < t) for the majority of spots in the sample, or that r/ro does not 
change very much in the same time interval so that the mean decay rate 
may be considered nearly constant in that period of time. 

For a time average taken over the whole decay phase, t± should be the 
time of area maximum, and £2 the end of the life of the spot. The number of 
spots in our sample for which both their birth and death were observed is, 
however, insufficient for a precise empirical determination of the distribution 
of log.D. 

One way to increase the number of spots in the sample is to define £1 and 
£2 as simply the first and last observations of a given spot, without respect 
to whether these coincide with the birth/death of the spot. However, we 
have just seen that for the departure from a Gaussian distribution to be 
significant it is important that we have £2 — £1 > r where r is the correlation 
time of decay rate fluctuations and that t/tq take a sufficiently wide range 
of values during the observation period. As these conditions are only valid 
for a relatively low fraction of all spots, the overall distribution derived in 
this manner may be expected to stay quite close to Gaussian. This may 
explain the findings of MMV93, where the method applied to determine D 
was indeed crudely the same as outlined in this paragraph. 

We therefore conclude that the probability distribution of (logD)' is 
close to normal, and we propose that it is this distribution that should be 
considered as fundamental both in a physical and in a mathematical sense. 
The approximately (but not exactly) Gaussian character of the distributions 
of log-D and log-D is then a secondary consequence of that fundamental 
distribution. 



3. The memories of a sunspot 

Figure 3 presents the autocorrelation of (log-D)' for sets of data pairs of 
different mean time separation (data pairs separated by 1,2,3,. . . consecutive 
observations) . It is apparent that the autocorrelation time (where the auto- 
correlation sinks below 1/e) is crudely ~ 3 , to be compared with a typical 
spot lifetime of ~ 7 in our sample. (This latter typical lifetime is entirely 
due to selection effects, but this is irrelevant for the present argument.) 

The existence of random time-dependent fluctuations in the decay rate 
with a correlation time shorter than the typical lifetime of spots leads to an 
apparent conceptual problem. As already mentioned in the Introduction, 
the mean decay law (1) has a curious property: it implies that sunspots 
have some "long-term memory" of their maximal radius even towards the 
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Figure 3. Autocorrelation of (logD)' for sets of data pairs of different mean time 
separation. 

end of their lives. Two small sunspots of identical size will have different 
expected daily decay rates if their maximal sizes were different. But if the 
decay process is subject to a random "noise" with a short correlation time, 
how is it possible for it to still possess a long-term "memory" in a statistical 
sense? 

The turbulent erosion model of sunspot decay (Petrovay and Moreno- 
Insertis, 1997, hereafter PM97) successfully predicted the law (1). In this 
model a "memory" effect is present due to the magnetic flux lost from 
the flux tube and piled up in its surroundings. The amount of this flux 
will influence the steepness of the field gradient just outside the current 
sheet (CS), which in turn determines the propagation speed of the CS 
(cf. Equation (10) in PM97). Thus, in the framework of this model any 
fluctuations in the decay rate must be due to fluctuations in the amount 
of the flux outside the tube. This flux can be considered the sum of two 
contributions: 

$ = $lost + $plage (6) 

where $i os t is the flux lost from the spot, while $ p i a ge is the flux of indepen- 
dent origin. This latter contribution is partly due to plage fields «£ p i,o which 
are originally present in the area, and partly to plage fields <3?' j which are 
carried into/out of the neighbourhood of the flux tube by horizontal flows, 
flux loss from other spots, and so on: 

$plage = $pl,0 + $pl (7) 

Both terms in Equation (6) show random fluctuations in time. Now av- 
eraging the equation we obviously have a relation for the expected value 
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of the flux around the spot (and therefore for the expected value of the 
decay rate) as a functional of the mean decay law (first term on the r.h.s. 
of Equation (6)) and as a function of the plage flux 3> p i o originally present 
in the site of the spot (second term; the other contributions <£' j to this 
term can be supposed to have zero mean). Specifically, if <3?pi ag e = 0, we 
are back to the models presented in PM97; the more general case will be 
treated in Section 5 below. The initial conditions thus fully determine the 
expected value of D at any later time. It is thus possible to have a long- 
term "memory" in the statistical sense while short-term fluctuations are still 
present. 

Clearly, the above reasoning also identifies a main physical factor con- 
tributing to the fluctuations in D: the ever-changing plage fields in the 
neighbourhood of the spot. Indeed, within the framework of the erosion 
model this is the only possible reason for deviations from the mean law. 
In reality, other contributions may originate in departures from the condi- 
tions of the erosion model (i.e. from depth- independence and axisymmetry) . 
However, as the erosion model proved to explain the mean decay law quite 
satisfactorily, it is not unreasonable to expect that the environmental effect 
of plage fields is at least one important factor in the origin of deviations 
from the mean decay law. 



4. Decay rate fluctuations and the plage field 

In the previous section we argued that one of the main physical factors 
contributing to the deviations from the mean decay law (1) is the envi- 
ronmental effect of plage fields in the neighbourhood of the spot. If this is 
the case then one should expect some correlation between the daily decay 
rate fluctuations and the plage flux surrounding the spot. In order to check 
this prediction, we determined the mean density B p \ of magnetic flux in 
an annulus of inner and outer radii r, and r around the spot from Kitt 
Peak magnetograms for each spot observation in our DPR data set (when 
a magnetogram made on the same day was available). 

As the magnetograms yield the line-of-sight field component only, this 
was deprojected into a vertical field by dividing it by cos a; a is the helio- 
centric angular distance from disk center. This method clearly introduces 
some extra scatter in the data but it cannot lead to systematic errors. The 
polarity is conventionally taken to be "positive" if identical to that of the 
spot in question. The results show little sensitivity to the precise values of 
n and r chosen; here we present results calculated with n = 1.5 ro and 
r = 3r . 

Figure 4a confirms our prediction of a correlation between B p \ and 
(logD)'. Indeed, the correlation coefficient is seen to differ from zero at 
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Figure 4- Scatterplots of the mean plage field strength B p i, in Gauss, around a spot vs. 
(logD)' and Ao, in MSH, for all measurements for spots that died on the visible hemi- 
sphere, excluding strongly non-axisymmetric spots (a, b, c), and for a typical simulated 
control data set, oi = a? — 0, 03 = 0.3 (d). Each panel is labeled with the value of the 
corresponding correlation coefficient. The fiducial line in (a) corresponds to the dashed 
line in Fig. 6 if B e — 400 G. The dashed and dotted lines in (c) correspond to spots with 
(logD)' =const. and lifetimes of 1 and 2 days, respectively. 



a confidence level well exceeding 2a. 1 Note that in this plot we have only 



1 The standard deviation of the correlation coefficient was here determined with the 
standard formula a — N^ 1 ' 2 , valid in the limit iV — > oo, so its value should be considered 
approximate. The alternative possibility of determining the significance level by a t- 
statistic was discarded as this method is only valid if the distribution of the data is 
binormal, possibly leading to gross errors otherwise. Given our ignorance regarding the 
origin and distribution of residual scatter in the B p j-(logD)' relation, the use of the 
approximate but robust N^ 1 ' 2 method seems safer, especially in the light of the recent 
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plage 




spots 

Figure 5. Sketch illustrating the origin of the correlation in Figure 3b. The annulus around 
small spots lies in the plage in its full extent, so these spots are surrounded by the stronger 
plage fields. 



included spots where the surrounding plage field was not too far from ax- 
isymmetric. The numerical criterion here was that, splitting the annulus in 

9 1 / 2 

five sectors, and calculating the mean field B{ in each sector, Bf [Bi < 1.2, 
where the bar indicates averaging over sectors. This selection actually only 
led to a slight improvement in the correlation which was quite similar for 
the complete sample. 

A doubt arises, however, concerning the origin of this correlation. It is 
apparent from Figure 4b that an even stronger correlation exists between 
B p i and Aq. This correlation can be explained by a trivial effect, illustrated 
in Figure 5: as the radius of the annulus is proportional to ro, for larger 
spots a large part of the annulus will lie outside the plage where mean 
fields are weaker, while for smaller spots the whole annulus will fall into 
the strong plage. On the other hand, owing to the observational selection 
effect treated in Section 3.3.2 of Paper I (spots with higher decay rates live 
shorter, and thus will be observed on fewer occasions), another correlation 
appears between log^o and (logD)' (Fig. 4c). The coupled effect of these 
two latter correlations may then be expected to give rise to a more indirect 
third correlation qualitatively similar to that seen in panel a. Thus the 
question arises if the anticorrelation seen between the plage field and the 
decay rate fluctuations is just the reflection of an observational bias, or it 
has some physical content? 

To resolve this dilemma, we generated artificial data sets with the same 
values of Aq and B p \ as the real data, but with (log-D)' values generated 
artificially by the formula 

(log D)' = a\ - a 2 B p i + a 3 x (8) 



controversion concerning the significance of correlations between time serii in the context 
of the solar neutrino problem (Oakley et al., 1994; Walther, 1998). 
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where x is a random variable of zero mean with a a = 1 Gaussian probability 
density. In order to test our null hypothesis we first set 02 = 0, i.e. assume 
that no real physical correlation exists between (logD)' and B p \. Then we 
submit our data set to a selection designed to mimic the observational bias 
affecting the real data. More specifically, we take the expected number of 
data for a given spot equal to its lifetime (as observations are normally made 
daily), which is in turn calculated assuming that (logD)' stays constant 
during the decay of the spot. 

The parameters a\ and 03 are freely adjustable, but subject to the con- 
dition that the resulting range of "observed" (log D)' values should coincide 
with the real range. After having computed a large number of simultaneous 
distributions with a wide variety of a± and 03 values we had to conclude that 
no parameter combination is able to yield a correlation coefficient higher 
than about 0.1 in typical realizations. The optimal agreement was found in 
the case a\ = ai = 0, 03 = 0.3, shown in Figure 4d. It is apparent that the 
correlation coefficient does not differ from zero by more than la. It seems 
therefore that the indirect effect of the correlations shown in panels b and c 
is insufficient to fully explain the correlation found in panel a, part of which 
must be a real physical effect. Indeed, if in our simulated data we allow 
02 7^ it becomes straightforward to find a good representation of the data: 
a\ = 0, ai = 0.3, aj, = 0.4 is for instance just such a case (not shown here). 

We thus conclude that while selection effects probably contribute to the 
anticorrelation of (logD)' with the plage field strength around spots, a part 
of this effect is likely to have a physical origin. 



5. Theoretical interpretation 

In order to understand the origin of a physical anticorrelation between 
(log D)' and the value of the external magnetic field strength here we extend 
the turbulent erosion model of sunspot decay for the case when the flux tube 
is embedded in a preexisting parallel homogeneous magnetic field. Apart 
from this generalization, we consider the same cylindrically symmetric case 
as in PM97, i.e. we solve the nonlinear diffusion equation 



M { " B) = I 



-Wf 



(9) 



where B is the magnetic flux density, f is the radial coordinate, and v is 
the turbulent magnetic diffusivity. For the dependence of the diffusivity on 
the magnetic field we still consider the same function satisfying the basic 
physical requirements: 

viB) = TTpkh' (10) 
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Figure 6. Current sheet velocity w as a function of the background field strength boo. O: 
a„ = 2.3, b = 30; □: a„ = 2.5, 6 = 10; A: a„ = 3.0, b = 10; *: a„ = 4.5, fe = 10; +: 
ct v — 7, bo = 3. The dashed line is just a fiducial straight line, drawn to aid the eye. Solid 
line: eq. (B-3) 



Here vq is the unperturbed value of the diffusivity, and a u is a parameter 
quantifying the steepness of the diffusivity cutoff near B e , the latter being 
the field strength where the diffusivity is reduced by 50%. Physically, one 
expects B e ~ B eq , the kinetic energy density of photospheric turbulence. 
Using the notation b = B/B e , our initial conditions will now be 



b(r) 



boo + 



bo ~ b c 



1 + (r/r ) 



22 



(11) 



where bo is the normalized field strength at the center of the tube, and boo 
is the external ("plage") field strength: we now assume 6oo 7^ 0. In this 
formula tq is the initial radius of the tube at time t = 0; for a comparison 
with observations it may be identified with the maximal equivalent radius 
of the spot. 

The values of the current sheet velocity w are plotted against boo for 
different numerical solutions of Equations (9)-(ll) in Figure 6. It is apparent 
that for boo = all the solutions are in agreement with the formula 



w ~ Cv /r (b - 1), 



(12) 



with C ~ 1.2, as expected. 



decay2pp.tex; 1/02/2008; 10:30; p. 12 



Making Sense of Sunspot Decay, II 13 

With intermediate values of a u ~ 2-3 the dependence of w on 600 is 
approximately logarithmic for 600 <; — 1, while for oi v > 4 it is strongly 
reduced for b^ > 0. The high-a^ behaviour may be interpreted in terms of 
a generalization of the dimensional argument leading to (12), see Appendix 
B. 

The actual value of a v is poorly known, but from Figure 6 it is clear that 
in the range — 1 < 600 < 1 the rate of flux loss should depend approximately 
logarithmically on b^, independently of the value of a v . In other words, the 
background field acts as multiplicative noise in the decay process. Setting 
B c ~ B cq ~ 400G, the photospheric value, this logarithmic relation is plot- 
ted over the observational data in Figure 4a. The apparent good agreement 
lends further support to the physical reality of the correlation found in the 
previous section. 



6. Conclusion 

6.1. Proposition: decay rate deviations as effects of the plage 

FIELD 

The theoretical results presented in Section 5 suggest the following scheme 
for the origin of random deviations from the mean parabolic decay law (1). 
It is rather natural to assume that the value of the background plage field 
strength of independent origin shows a Gaussian probability distribution 
around some expected value around a spot. Figure 6 then implies that the 
flux loss rate from the flux tube should be characterized by a lognormal 
distribution (just as we found in Section 2), and that an anticorrelation 
should exist between the logarithmic decay rate fluctuations (log-D)' and 
the strength of the plage field surrounding the spots (just as we found 
in Section 4). The generalization of the erosion model to the case of a 
tube embedded in an external homogeneous magnetic field then apparently 
offers a simple and elegant explanation for the presence of large lognormally 
distributed deviations from the mean decay rate in individual spots. 

It should be noted that a relation between the decay rate of sunspots and 
their environment was also suggested by Antalova and Macura (1985). While 
Moreno-Insertis and Vazquez (1988) correctly pointed out that the existence 
of a bimodal distribution of decay rates (i.e. well-separated branches of 
fast- and slow-decaying spots in the area-decay rate plane) in that study 
was the consequence of a trivial selection effect, the apparent physical 
difference found by Antalova and Macura (1985) between the two classes 
(slow-decaying spots tending to be "naked" spots in the sense defined by 
Liggett and Zirin, 1983) still remains, and it may possibly be related to the 
environmental effect proposed in the present work. 
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6.2. Caveats and call for more data 

However attractive the scheme outlined in the previous subsection may 
seem, it is almost certainly an oversimplification of the real situation. Firstly, 
it is obvious from Figure 4a that a large part of the scatter in (log D)' is 
independent of B p \. Though part of this scatter may be due to our method 
of deriving the vertical field component by deprojecting the line-of-sight 
component measured, etc., one may expect that effects like deviations from 
axial symmetry and sub-photospheric forces should necessarily contribute a 
great deal to deviations from the relation predicted by the erosional model. 
It is therefore certainly not the case that the lognormal distribution of decay 
rates can be explained by the plage effect alone: however, owing to the 
central value theorem of probability theory, the lognormal distribution can 
be saved if the other influences also act as multiplicative noise. 

Another caveat is related to the strong influence of the bias treated in 
Section 3.3.2 of Paper I throughout this analysis. Though, as we have seen, 
it is not impossible to correct for these effects, some fine points and simpli- 
fications in the correction process may always leave some doubt concerning 
the robustness of the conclusions derived. The main difficulties here are all 
due to the fact that our data (like all solar patrols realized in this outgoing 
century) were taken on a daily basis. One of the most emphatic conclusions 
of this study is that observations made on a daily basis are insufficient 
for an in-depth study of the sunspot decay process. What is needed is at 
least one observational campaign lasting several months, preferably at the 
time of solar maximum, with several simultaneous white-light heliogram and 
magnetogram observations a day. Only from the use of such higher quality 
data can we expect the final illumination of the sunspot decay problem. 
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Appendix A: A comparison of DPR and Mt Wilson sunspot 

areas 

In an attempt to extend our database we examined the possibility of 
incorporating Mt Wilson observations in our data (Howard, Gilman, and 
Gilman, 1984). One problem to be solved here was that in the Mt Wilson 
catalogue the spots are not identified from one day to the next, so in order 
to calculate decay rates the identification had to be done. We developed an 
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Figure 7. Total spot areas from the DPR 1977 plotted against their difference from the 
areas calculated by multiplying by 6.2 the corresponding umbral areas of the Mt Wilson 
catalogue. The size of the symbols is related to the central meridian (CM) distance of the 
spot, \L — Lcm\; the scale of symbol sizes is set by three examples given in the top right 
corner. 



algorithm for the identification of spots taking into account the differential 
rotation. A comparison of the results with the published DPR data, where 
available, showed that the identification has an acceptable rate of success, 
at least for the larger spots in each active region (which tend to dominate 
our data set). 

A further difficulty is caused however by the fact that the Mt Wilson 
data only yield estimated umbral areas. Umbral areas are notoriously more 
difficult to determine than umbra+penumbra areas (Gyori, 1998); besides, 
the method of determination applied was rather approximate (Howard, 
Gilman, and Gilman, 1984). While it is possible to translate umbral areas 
Ajj into total areas A by the simple linear rule A = 6.2 Ay (the value of the 
coefficient was determined from an analysis of DPR area data), this linear 
relation involves a large scatter, further contributing to the uncertainty of 
the area values. The resulting large errors are borne out in Figure 7, where 
we plotted the difference between the MtWilson area values thus derived 
and the corresponding DPR area values against DPR areas for 1977. It is 
apparent that the scatter is very large indeed, effectively excluding the use 
of Mt Wilson data for our present purpose. 

It may be worth noting that Figure 7 also seems to suggest a systematic 
trend: the areas computed from Mt Wilson data are lower than their DPR 
counterparts for large spots, while the reverse tends to be the case with small 
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Figure 8. Sketch showing the characteristic shape and scales of the magnetic field profile 
(solid) and of the stationary solution (B-l) (dashed) near the outer edge of the current 
sheet. r\ = b/(db/dr) at r = 1 



spots. One possible interpretation could be that the fixed ratio Aj A\j = 6.2 
used here may be incorrect, and that the umbral/total area ratio may in fact 
have some mean area dependence. However, a similar plot constructed using 
the DPR umbral areas instead of the Mt Wilson areas does not show any 
obvious trend (while it still shows a high scatter), which seems to exclude 
this explanation. The trend mentioned is then probably due to differences 
in the applied methods of area measurement. 



Appendix B: Analytical interpretation of the high-a t 
of the generalized erosion models 



behaviour 



For simplicity we set the units of time and length so that uq = 1 and 
6(1) = 1. Restricting our interest to the case &o ^ 1> from (12) w <C 1 
follows, so that in the range 1 < r < 1+n the solution may be approximated 
as the stationary solution of the linear diffusion equation in cylindrical 
geometry (cf. Fig. 8): 

lnr 



1 



(B-l) 



n 



From the numerical solutions we know that the erosional (i.e. nearly form- 
invariant propagating) solution forms after an initial transition time ti ~ 1. 
If the initial profile at t = is sufficiently steep [as is the case with (11)] 
then at t = U the amount of magnetic flux outside r = 1 equals the flux lost 
from within the CS, and the original plage flux there: 



2tt 



l+ri 



rb(r) dr ~ ivbo(r ~ 1) + ^00(^1 + 2 n) 



(B-2) 
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As we are still within about one diffusive timescale from t = 0, practically 
all this flux must still lie inside the radius 1 + r± , where the approximation 
(B-l) may be used for b(r). On the other hand, vq = 1 + wti, while w and 
r\ are related by w ~ [n(6o ~~ I)] -1 (PM97, eq. (10)), so, neglecting a term 
of order b^ 2 , we have 

(r? + 2n) (l + J- - bj) - ii±Ill! i n (i + n ) _ J^ n = (B-3) 

Setting ti = 1 and 6o = 10, Equation (B-3) can be solved numerically for n, 
and thus for w ~ [ri(&o — l)] -1 ; the solution turns out to be a very nearly 
linear function of boo (solid curve in Figure 6). 
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